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Abstract 

We study the asymptotic behavior of lattice power variations of two-parameter ambit 
fields that are driven by white noise. Our first result is a law of large numbers for such 
power variations. Under a constraint on the memory of the ambit field, normalized power 
variations are shown to converge to certain integral functionals of the volatility field asso- 
ciated to the ambit field, when the lattice spacing tends to zero. This law of large numbers 
holds also for thinned power variations that are computed by only including increments 
that are separated by gaps with a particular asymptotic behavior. Our second result is a 
related stable central limit theorem for thinned power variations. Additionally, we provide 
concrete examples of ambit fields that satisfy the assumptions of our limit theorems. 

Keywords: ambit field, power variation, law of large numbers, central limit theorem, chaos 
decomposition 
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1. Introduction 

1.1. Ambit fields and volatility 

A characteristic feature of many real-world random phenomena is that the magnitude or the 
intensity of realized fluctuations varies in time or space, or both. There are various terms 
used in different contexts that roughly correspond to this characteristic. To highlight two of 
them, in studies of turbulence, this is called intermittency, whereas in finance and economics 
the corresponding notion is (stochastic) volatility. Sudden extreme fluctuations — say, rapid 
changes in wind velocity or prices of financial securities — have often dire consequences, so 
understanding their statistical properties is clearly of key importance. 

Barndorff- Nielsen and Schmiegel [11, 12] have introduced a class of Levy-based random 
fields, for which they coined the name ambit field, to model space-time random phenomena that 
exhibit intermittency or stochastic volatility. The primary application of ambit fields has been 
phenomenological modeling of turbulent velocity fields. Additionally, Barndorff-Nielsen, Benth, 
and Veraart [2] have recently applied ambit fields to modeling of the term structure of forward 
prices of electricity. Electricity prices, in particular, are prone to rapid changes and spikes since 
the supply of electricity is inherently inelastic and electricity cannot be stored efficiently. It 
is also worth mentioning that, at a more theoretical level, some ambit fields have been found 
to arise as solutions to certain stochastic partial differential equations [3]. Barndorff-Nielsen, 
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Benth, and Veraart [4] provide a survey on recent results on ambit fields and related ambit 
processes. 

In this paper, we study the asymptotic behavior of power variations of a two-parameter 
ambit field driven by white noise, with a view towards measuring the realized volatility of the 
ambit field. Specifically, we consider ambit field (X(s.t))o,t)e[o,i] 2 ; defined via the equation 

Yu t )= g{s-u,t~v)a {U:V) W{du,Av), (1.1) 

J A(s,t) 

where the integrator W is a white noise on R 2 and the integrand is defined in terms of a positive- 
valued, continuous volatility field (c( s ,i))(s,t)eR 2 an d a weight function g G L 2 (M. 2 ). The integral 
in (1.1) is computed over the set A(s,t) C M 2 , which is known as the ambit set associated to 
the point (s,t). More figuratively, A(s, t) defines the "ambit" of noise and volatility innovations 
that influence Yi s ^y We use here the common specification of A(s,t) as a translation of some 
fixed Borel set A C R 2 , viz., 

A(s,t) :=A+(s,t) := {(u + s,v + t) : (u,v) G A}. (1.2) 

The shape of the set A has a strong influence on the probabilistic properties of Y. When the 
parameter t is interpreted as time, it is customary to assume that A C K x (— oo,0], so that 
only past innovations can influence the present. We refer to [3] for a discussion on the possible 
shapes of A in various modeling contexts. We consider here only the case where the volatility 
field a and the white noise W are independent. In this case the integral in (1.1) can be defined in 
a straightforward manner as a Wiener integral, conditional on a. (Ambit fields with volatilities 
that do depend on the driving white noise can be defined, but then the integration theory 
becomes more involved, see [3] for details. Moreover, the general framework of ambit fields also 
accommodates non-Gaussian random measures, Levy bases, as driving noise.) 

The power variations we study are defined over observations of Y on a square lattice in 
[0, l] 2 using rectangular increments (see Section 2.3 for precise definitions). The spacing of the 
square lattice is 1/n, and we let n — > oo in the asymptotic results. In addition to ordinary 
power variations that involve all of the available increments, we consider also thinned power 
variations that are computed using only every k n -th increment in the lattice. Asymptotically, 
we let k n — > oo so that k n /n —> 0. Similar procedures have been considered in the context 
of Gaussian processes by Lang and Roueff [25] and, more recently, in the context of Brownian 
semistationary processes by Corcuera et al. [19]. 

Our first result is a functional law of large numbers for both ordinary and thinned power 
variations (Theorem 2.6). Under an assumption that constrains the memory of Y through the 
so-called concentration measures associated to the weight function g (Assumption 2.5), we show 
that the suitably scaled power variation of Y converges in probability to an integral functional 
of the volatility field a. Under a more restrictive and quantitative version of Assumption 2.5 
(which appears as Assumption 2.8), we also obtain a stable functional central limit theorem 
for thinned power variations (Theorem 2.12) with a conditionally Gaussian random field as the 
limit. We give some explicit examples of weight functions g that satisfy Assumptions 2.5 or 2.8 
in Section 2.5. 

The motivation of this paper is twofold. On the one hand, the study of the asymptotics 
of power variations of ambit fields is interesting from a probabilistic perspective, as it provides 
information on the fine structure of the realizations of ambit fields. On the other hand, in 
practical situations it is of interest to draw inference of volatility statistics of the form 
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for p > 0, based on discrete observations of the ambit field Y. Our law of large numbers 
establishes a sufficient condition that the suitably scaled p-th power variation of Y over [0, s] x 
[0, t] converges to (1.3). This could be seen as a first step towards a theory of volatility estimation 
for ambit fields. 

1.2. Related literature 

There is a wealth of literature on laws of large numbers and central limit theorems for power, 
bipower, and multipower variations of (one-parameter) stochastic processes. Notably, semi- 
martingales are well catered for, see the monograph by Jacod and Protter [22] for a recent 
survey of the results. Similar results for non-semimartingales are, for obvious reasons, more 
case-specific. Closely relevant to the present paper are the results for Gaussian processes with 
stationary increments [5, 8] and Brownian semistationary processes [6, 7, 19]. In fact, a Brow- 
nian semistationary process is the one-parameter counterpart of an ambit field driven by white 
noise. The proofs of the central limit theorems in [5, 6, 7, 8, 19] use method that involves 
Gaussian approximations of iterated Wiener integrals, due to Nualart and Peccati [29]. We 
employ a similar approach, adapted to the two-parameter setting, in the proof of our central 
limit theorem. 

Barndorff- Nielsen and Graversen [10] have recently obtained a law of large numbers for the 
quadratic variation of an ambit process driven by white noise in a space-time setting. The 
probabilistic setup they consider is identical to ours, but their quadratic variation is defined 
over observations along a line in two-dimensional spacetime, instead of a square lattice. The 
proof of our law of large numbers is inspired by the arguments used in [10]. 

Asymptotic results for lattice power variations of random fields with two or more parameters 
are scarser. Kawada [24] proves a law of large numbers for general variations of a class of Gaus- 
sian random fields, extending an earlier result of Berman [14]. A related central limit theorem, 
which covers only quadratic variations, is due to Deo [21]. More recently, Reveillac [32, 33] 
has obtained central limit theorems for weighted quadratic variations of ordinary and fractional 
Brownian sheets. Similar results, which include also non-central limit theorems, applying to 
more general Hermite variations of fractional Brownian sheets appear in the papers by Breton 
[18] and Reveillac, Stauch, and Tudor [34]. 

2. Definitions and main results 

2.1. Notation 

For any z G R 2 , non-empty A C M 2 , and r > 0, we write B(z,r) := {( e K 2 : \\C - z\\ < r}, 
and A r := U^gA -^(C r )- Moreover, A stands for the closure of A in M 2 . 

For any s, t € K., we use s A t := min(s, t) and s V t := max(s, t), [sj '■— max{r G Z : r ^ s}, 
\s~\ := min{r G Z : r ^ s}, and {x} := x — \_x\ . It will be convenient to write s <e t 
(resp. s >e t) whenever there exists Cg > that depends only on the parameter 9, such that 
s ^ Cgt (resp. Cgs ^ t). We write s x# t to signify that both s <g t and s >g t hold. 

We denote the weak convergence of probability measures by A, the convergence of random 
elements in law by A, and the space of Borel probability measures on M 2 by ^(IR 2 ). The 
support of v G "P(1R 2 ), or briefly suppi^, is the smallest closed set with full ^-measure, given by 

f) r >o{ z e R2 : v{ B ( z , r )) > °}- The Lebes 

gue measure on M. d is denoted by and the Dirac 

measure at z G M. d by S z . 

For any q > 0, we write m q := E[|X| 9 ], where X ~ N(0, 1). Finally, \A\ stands for the 
number elements in a finite set A, and we use the conventions N := {1, 2, . . .} and No := NU{0}. 
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2.2. Rigorous definition of the ambit field 

Let W be a white noise on [— 1, l] 2 with the Lebesgue measure A2 as the control measure. Recall 
that this means that W is a zero-mean Gaussian process indexed by B{[— 1, l] 2 ) with covariance 
F,[W(A)W(B)] = X 2 (A n B) for any A, B e 8([-l,l] 2 ). Throughout this paper, we consider 
an ambit field Y given by 

y (s , t) := J g(s-u,t- v)a {UjV) W(du, dv), (s, t) G [0, l] 2 , (2.1) 

where g G L 2 (R 2 ) is a non-vanishing weight function and (o7 S) t))(s,t)e[-i,i] 2 i s a continuous, 
strictly positive volatility field, independent of W. Let us denote by A g the essential support of 
g (see, e.g., [26, p. 13] for the definition). In (2.1) it suffices to integrate over the set 

A(s,t) := -A g + (s,t). 

Thus, we recover the setting outlined in (1.1) and (1.2) with A — —A g . To ensure that A(s, t) C 
[-1, l] 2 for all (s,t) e [0, 1], we assume that A g C [0, l] 2 . 

The stochastic integral in (2.1) is a conditional Wiener integral with respect to W, defined 
as follows. Due to the independence of a and W, we may assume without loss of generality that 
the underlying probability space is the completion of the product space 

(fl w x T w ® F a ,P w ® Pa), 

where (fiwj-FwjPw) carries the white noise W, so that Tw — a{W{A) : A € B([—l, l] 2 )}, and 

(fi„,.F«r,P«,) - (C([-l,l] 2 ,M + ),£(C([-l,l] 2 ,]R + )),P ff ) 

is the canonical probability space of cr, i.e., ar xt \(u)) := uj(s,t) for any u) € fl a and (s,t) £ 
[— 1,1] 2 . Then, for any uj £ Q a and (s,i) e [— 1,1] 2 , we define Y( s t )(-, w) to be the Wiener 
integral of the function (u, v) g(s — u, t—v)ui(u, v), which belongs to L 2 ([— 1, l] 2 ), with respect 
to W. Since the Wiener integral is a linear isometry between the integrand space L 2 ([— 1, l] 2 ) 
and L 2 (Q,y/) (see, e.g., [28, pp. 7-8] for details), no issues will arise with the measurability of 

Remark 2.2. Given any continuous function L : [0, 1] — > [0, l] 2 , i.e., a curve, we may define a 
stochastic process (^r(t))te[o,i] > gi vm g the description of the ambit field Y as seen by an observer 
moving along the curve L. Such processes are called ambit processes. Barndorff-Niclsen and 
Graversen [10] study the limit behavior of the quadratic variation of (^r(t))*G[o,i] m t ne case 
where T is a line segment, establishing sufficient conditions for the law of large numbers. 

2.3. Power variation and concentration measure 

For a two-parameter random field, an increment is naturally defined over a rectangle in the 
parameter space. Specifically, the rectangular increment of the ambit field Y over R := (si, S2] x 
(ti,t 2 ] C [0, l] 2 is defined as 

Y(R) := Y {S2M) - Y (siM) - Y {s2ttl) + F (si , 4l) . (2.3) 

The definition (2.3) is standard in the literature of random fields, and can be recovered for 
example by partial differencing of Y( s t j with respect to s and t — or vice versa. Although 
not needed in the sequel, it is worthing pointing out the fact that the map R H> Y(R) can be 
extended to a finitely additive random measure on the algebra generated by finite unions and 
intersections of rectangles in [0, l] 2 , which motivates the notation Y(R). 
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For fixed p > 0, we shall consider the p-th power variation of Y over the square lattice 
S n := {(^, : i, j = 0, 1, . . . ,n} C [0, l] 2 for any n G N. Based on the values of Y on the 
lattice 5„, we may compute the increments of Y over the rectangles 

R tij) := ((* ~~ ^AM/"] x ((i ^ i)/ n >J'/ n L i, j* = i,...,n. 

Using them, we define the p-th power variation of Y" over S n by 

V$ } (k,n):= £ £ l y «l*i))r. (M)e[0,l] 2 . (2-4) 

»=1 3=1 

where fc e N is a thinning parameter. This allows us to take only every /c-th increment into 
account when computing the power variation. The case k — 1 corresponds to ordinary power 
variations whereas letting k > 1 gives rise to thinned power variations. Note that we regard 
n) as a random field on [0, 1] . 
To state the assumptions of our limit theorems, we need to introduce a technical device that 
controls the interdependence of the increments appearing in (2.4). Let us first define h n G L 2 (M. 2 ) 
for any n £ N by 

h n (s,t) := g(s,t) - g(s - l/n,t) - g(s,t- 1/n) + g(s - l/n,t- l/n), 
which, in fact, enables us to write succinctly 

Y(R^\- ) ) = f h n (i/n- u,j/n- v)cr( UiV )W(du,dv). 

Since g is non- vanishing, we have c„ := J R2 h n (z) 2 dz € (0, oo). Thus, we may define 7r„ e 7 , (M 2 ) 
by 

ir n (dz) := Tr n (z)dz, where ir n {z) := hn Z . 

The probability measure 7r„ is a so-called concentration measure, analogous to the ones appearing 
in earlier papers on ambit processes [10, p. 265] and Brownian semistationary processes [6, 
p. 1166]. Roughly speaking, the strength of the interdependence of the increments (2.4) is 
related to how dispersed ir n is. Our limit theorems are based on the key assumption that the 
interdependence is not "too strong" , in the sense that the sequence 7Ti , 7T2, . . . converges weakly 
to a probability measure that is supported on a "small" subset of M 2 . 

2.4. Limit theorems 

We state now the main results of the paper. Their proofs, along with some auxiliary lemmas, 
are deferred to Sections 3 and 4. In this section, (fc n )„gN stands for a non-decreasing sequence 
of natural numbers, which we shall use as the values of the thinning parameter, such that 
£n : = k n /n —> 0. However, the assumption that k n — ¥ oo is not imposed yet. 

Our first result is a functional law of large numbers for V^(k n ,n). The key assumption, 
which was alluded to above, behind the law of large numbers is the following. 

Assumption 2.5. There exists it e V(R 2 ) such that A2(supp7r) = and n n A n. 

The condition A2(supp7r) = holds, for example, when tt is concentrated on a curve and, 
in particular, when it is a convex combination of finitely many Dirac measures. Examples of 
weight functions g that satisfy Assumption 2.5 are given in Section 2.5. 

In the statements below, D([0, l] 2 ) C E' 0,1 ' stands for the natural two-parameter general- 
ization of the cadlag space D([0, 1]) C M^ 0,1 !. We endow this space with the uniform topology. 
Appendix B recalls the precise definition of D([0, l] 2 ), along with some useful related facts. 
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Theorem 2.6 (Law of large numbers). If Assumption 2.5 holds, then 
4^'(^ 1 ") P - mpE^'i C([0,1]) 2 , 

where 

^ := U af u _^_ T) n(d^r)Y dudv, (s,t) £ [0,1] 2 . 

Remark 2.7. Assumption 2.5 is slightly more restrictive than mere mutual singularity of 7r and 
A 2 . Indeed, the proof of Theorem 2.6 uses a separation argument that relies on the existence of 
a closed A2-null set with full 7r-measure. 

The case where ir = 5( S0) t ), for some (so,to) € [0, l] 2 , is of particular interest. Then, we 
have 

rs—s ft— to 

S p ; = / / a? ,dwh/. 

(s,t) / / (u,t>) 

•»— «o " / — *o 

From a practical point of view, the case where ir is not & Dirac measure is somewhat undesirable. 
Then the random field E^' 71 ^ "sees" merely a weighted space-time average of a, and inferring 
the "pure" a may become impossible. 

Our second result is a functional central limit theorem for V^(k n , n). Here, we concentrate 
on the case where 7r is a Dirac measure and k n — > oo. For the needs of the central limit theorem, 
we refine the Assumption 2.5 by quantifying the speed of the convergence 7r„ — » ir as follows. 

Assumption 2.8. There exist open sets Ei, E^, . . ■ cK 2 and Zq := (sq, to) £ [0, l] 2 such that 
for all n £ N, 

(i) Z £ E n , 

(ii) \ 2 {E n n (E n + (s, t))) = for any (s, t) £ R 2 such that |s| V |t| ^ £„, 

(iii) ^„(IR 2 \£;„) = o(e 2 ). 

The sets E\,Ei,... should be seen as shrinking "neighborhoods" of the point Zq. In fact, 
items (i) and (ii) imply that for all n £ N, 

E n C [s 

Thus, by item (iii), Assumption 2.5 holds with ir = 8 ZQ . 

The central limit theorem is stated in terms of stable convergence in law, a notion due to 
Renyi [31], which is the standard mode of convergence used in central limit theorems for power, 
bipower, and multipowcr variations of stochastic processes. For the convenience of the reader, 
we recall here the definition. 

Definition 2.9 (Stable convergence in law). Let U\,U2,... be random elements in a metric 
space U, defined on the probability space (fijJ 7 , P), and let U be a random element in U, 
defined on (f2', J 7 ', P'), an extension of (ft, J 7 , P). When Q C T is a c-algebra, we say that 

XJ\, U%, . .. converge Q -stably in law to U and write U n — \ U, if 

B[f(U n )V] > E'[f(U)V] (2.10) 

n— too 

for any bounded, ^/-measurable random variable V and bounded / £ C(U,R). 

Remark 2.11. Choosing V = 1 in (2.10) shows that stable convergence implies ordinary conver- 
gence in law. However, the converse is not true in general. 
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Theorem 2.12 (Central limit theorem). If Assumption 2.8 holds, then 



^(yW^.nJ-E^^'^.n)]) _^( m2p _ m 2)i/2 3 (p) mD^if), ( 2 .13) 

/ ^^(du.dw), (M)e[0,l] 2 



"(p) ._ 
-( S) t) — 

[-s ,s-s ] x [-t ,t-t ] 

and W 1 - is a white noise on [0, l] 2 with control measure X 2 , independent of J 7 , defined on an 
extension of (f2, J-, P) . 

2.5. Weight functions 

We shall now briefly discuss some examples of weight functions g that satisfy Assumptions 2.5 
or 2.8. 

2.5.1. Uniform weight function 

Perhaps the simplest possible weight function is such that it assigns uniform weight over a 
rectangle. More concretely, let 

.9 ; = l[si,S 2 ]x[tl,t 2 ]! 

where ^ Si < s 2 ^ 1 and ^ t\ < t 2 ^ 1. For any n > 1/ ((s 2 — s i) A (t 2 — ^i)), we have 

h n = l[si,si+l/n]x[ti,ti+l/n] ~ l[s 2 ,S2+l/n]x[ti,ti+l/n] 

— l[s 1 ,s 1 + l/n]x[t 2 ,t 2 + l/n] + l[s 2 ,s 2 + l/n]x[i 2 ,i 3 +l/n] 

almost everywhere. It is easy to check that then c„ = 4/n and that Assumption 2.5 holds with 
7r = (l/4)(<5( Slitl ) + 6( S2 . tl ) + S( Sut2 ) + 8( S2 ,t 2 ))- Thus, Theorem 2.6 implies that 

V(2)(1 ' n) ^^ / / ( a (u-s u v-H) + 4»-«».«-*0 + o-f u - suv - t2 ) +o-f u _ S2 ^ t2) )dudv 
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in -D([0, 1] ). Assumption 2.8, of course, cannot hold under this specification of g. 
2.5.2. Weight function with a singularity 

To satisfy Assumption 2.8, the weights imposed by g should be concentrated to a neighborhood 
of some point in [0, l] 2 . For example, let us consider g £ L 2 (M. 2 ) with a singularity at zero, given 

{(.vtr^vt), (M)e(0i) 2 , 

where a € (0, 1) and ^ G C x (0, 1) is such that lim s ^ + £( s ) + °! linis-n- ^(s) = 0, and ||f := 
sup sg ( 01 ) \£'(s)\ < 00. Note that, necessarily, we have also p||oo : = su P s e(o,i) l^( s )l < 00. A 
simple example of such a function is i(s) := 1 — s. 

Assumption 2.8 holds under this specification provided that the thinning parameter k n has 
suitably fast rate of growth. The following result gives a sufficient condition in terms of the 
asymptotic behavior of e n . Its proof is carried out in Section 5.1. 

Proposition 2.15 (Weight function with a singularity). Suppose that g is given by (2.14). 
(1) Assumption 2.5 holds with n = 5q, 
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(2) If e n x n~~ K , where < k < a when a € (0, 1/2) and < k < (2a + l)/(2a + 3) w/ierc 
a G [1/2, 1), i/ien, Assumption 2.8 holds with Zq — and £"„ = (0,e„) 2 . 

Remark 2.16. If we assume further that lim ;s _ i . 1 _ l'(s) ~ 0, then it is possible to show that 
< k < (2a + l)/(2a + 3) is a sufficient condition for all a € (0, 1). 

2.5.3. Weight function supported on a triangle 

As another example, let a G (1/2, 1) and £ as above, and define g G L 2 (M. 2 ) through 



t- a t(t), (s,t)eT, 
0, (s,t)eM. 2 \T, 



=\ n , ... 2lm ( 2 - 17 ) 



where T := {(s,i) : (1 — t)/2 < s < (1 + t)/2, < t < 1} is the isosceles triangle with vertices 
(1/2, 0), (0, 1), and (1, 1). Such a weight function is typical in space-time modeling of turbulence 
(see, e.g., [9, 11, 35]). Interpreting s as a one-dimensional space variable and t as time, the set 
T (or more appropriately —T) can be seen as a causality cone. 

Due to the different shape of the support, under this specification of g we require that the 
thinning parameter grows at a faster rate compared to the preceding example. The proof of the 
following result is very similar to the one of Proposition 2.15, so it is merely sketched in Section 
5.2. 

Proposition 2.18 (Weight function supported on a triangle). Suppose that g is given by (2.17). 

(1) Assumption 2.5 holds with ir = 5 Zo , where z = (1/2,0). 

(2) If £„ x n~ K , where k G (0, (2a — l)/(2a + 1)), then Assumption 2.8 holds with E n = 
(1/2 - e n /2, 1/2 + e B /2) x (0,e„/2). 

Remark 2.19. It is evident from the proof that Proposition 2.18 can be easily extended to a 
weight function g whose essential support is a "small perturbation" of the triangle T . 

2.6. Some comments on the results 

2.6.1. Measurement of relative volatility 

A practical difficulty in using Theorem 2.6 is that the power variations need to be scaled ap- 
propriately and the scaling depends on the unknown weight function g and may be difficult to 
compute precisely. In fact, it is evident that the volatility field a cannot even be determined un- 
ambigously unless g normalized a priori. However, often we are more interested in the variation 
a rather than its precise level, which may not, thus, be very informative due to the ambiguity 
caused by the lack of normalization. It is key to note that the variation of a is captured also by 
the relative integrated volatility field 

J .i f i { 2 \ - , (M)e[0,l] 2 . (2.20) 

Jo Jo a (u,v) dudv 

Quantities of the form (2.20) can be obtained as the limits of certain ratios of (unsealed) power 
variations, which are statistically feasible. More precisely, Theorem 2.6 readily implies that for 
any p > 0, 

V ( [%(l,n) ™ £ of^dvdv 
The use of relative volatility statistics, in general, is elaborated in a forthcoming paper [13]. 
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2.6.2. Bias in the central limit theorem 

Note that in Theorem 2.12, the scaled power variation e 2 c _p / 2 ^ p )(fc n , n) is centered around 
its expectation e 2 c~ p / 2 Ew \V.(k n , n)] , instead of the limit m p £( p ' 7r ) given by the law of large 
numbers. While it is shown in the proof of Theorem 2.6 that, under Assumption 2.5, 

2 

^E ff [Vg>)(fen, n)] m p sj^ for any (s, i) e [0, l] 2 , (2.21) 

the rate of convergence in (2.21) appears to be in most, if not all, cases too slow that we could 
replace s 2 n c' p ^'E w [V {p) (k n ,n)] with m p Y,^ in (2.13). 

An asymptotically non-negligible bias is present even in the most well-behaved case with 
constant <7( Sjt ) = o. Namely, we have then 

= - m ^({i} t+ {^} s+o(1) )- 

One can show that for almost any (s,t) £ [0, l] 2 , 

limsupfj — } + { — }) >0, 

and, consequently, 

2.6.3. Extending the central limit theorem beyond thinned power variations 

Is it possible to extend Theorem 2.12 to cover ordinary power variations? Quite possibly, but 
we expect that the limit would not remain the same. In fact, we conjecture that the situation is 
analogous to Brownian semistationary (BSS) processes (see [19]). Recall that ordinary power 
variations of BSS processes, under certain conditions, satisfy a central limit theorem [19, The- 
orem 3.2] with a limit analogous to S^ p \ but multiplied with a constant that is strictly larger 
than (rti2 V — m 2 ) 1 / 2 , whereas the limit in the corresponding result for thinned power variations 
[19, Theorem 4.4] has the factor {mi P — m 2 ) 1 / 2 . This is a consequence of the non-generate 
limiting correlation structure (which identical to the one of fractional Brownian noise) of the 
increments of a BSS process. Thinning decreases the asymptotic variance in the central limit 
theorem through "decorrelation" of the increments, but at the expense of rate of convergence. 

While our Theorem 2.12 is analogous to Theorem 4.4 of [19], obtaining a central limit theorem 
for unthinned power variations, akin to Theorem 3.2 of [19], is currently an open problem, which 
we hope to address in future work, along with allowing for a that depends on the driving noise. 
The key problem is the identification of the limiting correlation structure of the increments. 
However, it seems that such a result cannot be accomplished by a straightforward modification 
of the arguments in [6] since the one-dimensional regular variation techniques used with BSS 
processes appear unapplicable in our setting due to the additional dimension. We also expect 
that, like in [6, 7, 19], such a result would require stronger assumptions on the dependence 
structure of the ambit field — beyond what we formulate using the concentration measures — 
and a smoothness condition on a. 
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3. Law of large numbers 



In this section, we prove the law of large numbers for power variations, Theorem 2.6. The proof 
is based on the conditional Gaussianity of the ambit field Y and, in particular, on a covariance 
bound for nonlinear transformations of jointly Gaussian random variables, which we will review 
first. 

3.1. Hermite polynomials and a covariance bound 

Recall that the Hermite polynomials Hq, Hi, H 2 , ■ ■ ■ on R are uniquely defined through the 
generating function 

( t 2 \ °° 
exp [tx - — = V t n H n (x), ieE. 

V 2 / „=i 

They are orthogonal polynomials with respect to the Gaussian measure 7 on K. More precisely, 
if (Xi,X 2 ) is a Gaussian random vector such that E[X X ] = E[X 2 ] = and ~E[Xf] = ELYf] = 1, 
then (cf. [28, Lemma 1.1.1]) 

n\B[H n (X 1 )B m (X 2 )] = i 1 1 2J ' ' (3.1) 

10, nf m. 

Thus, \yr)}.H n : n £ No} is an orthonormal basis of L 2 (IR,7) and, in particular, for any 
/ e L 2 (R,j) there exists (a ,ai, ...)(£ ^ 2 (N ) such that 

00 

f = J2 a n^-H n inL 2 (M, 7 ). (3.2) 

n=0 

The index of the leading non-zero coefficient in the expansion (3.2), that is, min{fc S No : at ^ 
0}, is known as the Hermite rank of the function /. 

Using (3.1) and (3.2), it is straightforward to establish the following bound for covariances 
of functions of jointly Gaussian random variables. This simple inequality is, in fact, a special 
case of a far more general result due to Taqqu [37, Lemma 4.5]. 

Lemma 3.3 (Covariance). Let (Xi, X2) be as above. If f € L 2 (M.,-f) has Hermite rank r € N, 
then 

\E[f(X 1 )f(X 2 )}\ < f \E[XiX 2 ]\ 9 for any q e [0,r]. 

For any p > 0, write u p (x) := \x\ p — m p , x € M. Clearly, m p € L 2 (K,7) and Gaussian 
integration by parts shows that the Hermite rank of u p is 2. Thus, Lemma 3.3 implies that 

[CovOXrl^, |X 2 H| < p lE^X^I" for any g e [0,2], (3.4) 

which will be instrumental in the proof of Theorem 2.6, below. 

3.2. Proof of Theorem 2.6 

Prior to proving Theorem 2.6, we still need to establish a simple fact that follows from the 
convergence 7r„ tt. To this end, recall that the Levy-Prohorov distance of /1, v € ^(K 2 ) is 
defined as 

d(At,i>) := inf {e > : fx(E) < ^(^ e ) + e, i/(£T) < /i(75 e ) + e for all E e /3(M 2 )}. 

The Levy-Prohorov distance is a metric on ^(K 2 ) and 7r„ 7r holds if and only if d{ir n , tt) — > 
(see, e.g., [16, p. 72]). Below, we write i? := supp7r, for the sake of brevity. 
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Lemma 3.5 (Concentration). If 7r n n, then there exists positive numbers (a n ) such that 
a n I and ir n (B an ) — > 1. 

Proof. Let (a„) be such that a„ J, and a n > d(7r„,7r) for any n G N. By the definition of 
the Levy-Prohorov distance, tt(B) ^ ir n (B a ™) + a n for any n £ N. Since ir(B) — 1, we have 
ir n {B a ") > l-a„ -> 1. □ 

Proof of Theorem 2.6. Clearly, we have V^ p ^(fc„,n) V^^(fc n ,n) if s it and t ^ v. Thus, 
by Lemma A.l, it suffices to establish pointwise convergence 

r s r t / r \ p/2 



2 -P/2 T/ (P) 



V (% {k ^ n) m P [I (J ''l-^-^^dr^ dudv (3.6) 

for any (s,i) G [0, l] 2 . More precisely, we show (3.6) conditional on the realization of a. Under 
this conditioning, we may regard a as a non-random element of C([— 1, l] 2 , M+) and V as a 
Gaussian random field. 
Let us first show that 



lhn elc-^ w [V^ t) {k n ,n)] = m p J J* (J <jf u _^_ T) n(d^ drjY' dudv. 



(3.7) 



Since 



TOpC " /2 (/ <7 ( i /™-«.i/™- T ) 7rn ^' dr )) ' 



we have 

Ls/e„J |t/e„J / „ x p/2 



i=l 3=1 ^ J 

(/ tfWn-^Hn-T^nidt^T)^ dudV, 



o Jo 



where \x\ n :— e n \x/e n ~\ and [x\ n := e n \x/e n \ for any x G R and n G N. Since [sjn — > s 
and [t\ n — > t as n — > oo, the convergence (3.7) follows from Lebesgue's dominated convergence 
theorem, provided that for any (u, v) G [0, l] 2 , 



ii^L J ^CMnAi-a^UAi-T^nCd^dr) = / cr ( „_ ? „_ r) 7r(d£, dr), 

which, in turn, is a straightforward consequence of the uniform continuity of the realization of 
o~ and the convergence ir n —> n. 

Now, (3.6) follows from Chebyshev's inequality, provided that 

lim 4c"f Var w [V ( % (k n , n)] = 0. (3.8) 

To show (3.8), we expand 
e^Var w [V ( J] ) (fc n ,n)] 

= 4 E J E J c -*cov w [|y .^IMy^^n. (3-9) 

*i,*a=l Jij2=l 
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Using the inequality (3.4) and the relation 

E Wl y KS))| 2 ] ~- c »> »,j = l,...,n, neN, (3.10) 



we obtain 



C ;1Cov w [|y(i?g iiifcnii) )|^|y(i?£ i2ifeni2) ) 



<^c„|E w [F(i?g ii!fcnii) )y(iij£ i2ifcn . 2) )] 



7i"n(^, t) 1/2 tt„(£ + e„(ii - i 2 ), r + e„(ji - j 2 )) 1/2 d£di 
Applying this bound to (3.9), we arrive at 



s t s t 



£ n C r/ Var w[V( ( s P i)(^,™)] ^cr,p j j j J ^n{ui,Vi,U2,V 2 )du 1 dv 1 du 2 dV2, 



where 

U n (u 1 ,v 1 ,u 2 ,v 2 ) ■■= / ^ n {i,T) 1/2 'k n {i+ \u{\ n - \u 2 ~\ n ,T + \v{] n - [w 2 ln) 1/2 dfd7 



The Cauchy-Schwarz inequality ensures that the functions IIi,Il2, • • ■ are uniformly bounded 
on ([0, s] x [0, t]) 2 . Thus, by Lebesgue's dominated convergence theorem, it suffices to show that 
n„ tends to zero almost everywhere as n — > oo. We will split this task into two parts by treating 
separately 

Ii^\ui,vi,u 2 ,v 2 ) := / 7r„(£,T) 1/2 7r„(£ + \u{\ n - \u 2 ] n ,T+ \v{]„ - \v 2 ] n) 1/2 d£dr 

and 

U^(ui,vi,u 2 ,v 2 ) := / 7r n (^,T) 1/2 7r n (^ + \u{\ n ~ \u2]n,r + \v{] n - \v 2 ] n ) 1/2 d£dr, 

where (a n ) is a sequence of positive real numbers such that a n J, and ir n (B an ) — > 1, the 
existence of which is ensured by Lemma 3.5. Applying the Cauchy-Schwarz inequality to IT 1 ), 
we obtain 

Ui 1 '>{u 1 ,v l7 u 2 ,v 2 ) 2 < n n (R 2 \B a ~) [ 7r„(£ + \u x ] n - \u 2 ] n ,T + \v{\ n - \v 2 ] n )d£dT 

SC 1-7T„(B Q ") ► (). 

n— yoo 

Similarly, in the case of n( 2 ) we obtain 

I^ 2 \ui,vi,u 2 ,v 2 ) 2 ^ / 7T„(£ + \ui\ n - \u 2 ~\ n ,T + \v{] n - [v 2 ~|„)d£dT, (3.11) 

where, however, a slightly more elaborate argument, inspired by the proof of Lemma 1 in [10], 
is needed to show convergence to zero. 

By Urysohn's lemma, for any 5 > there exists ip$ £ C(R 2 , [0, 1]) such that tpg(z) = 1 for 
z £ W and <p s (z) = for z £ R 2 \ B 2S . From (3.11) we deduce, thus, 

limsupr4 2) (ui,ui,u 2 ,u 2 ) 2 lim / <p s (£ , r)7r„(£ + \u{\„ - \u 2 ] n ,r+ \vi] n - \v 2 ~\ n )d£,dT 

n— >oo n— >-oo J 

= lim / tp s (£+ \u2~\n - \ui~\n,T+ \v 2 ] n - [vi~| „W n (d£, d-r) 

71— ¥00 J 

tps(£, +U 2 -U!,T + v 2 - wi)7r(d^,dr), 
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where we used the bound | \x] n — x\ < e n , for all x £ R and n £ N, and the observation that ip E 
is, in fact, uniformly continuous. Since tp$ converges pointwise to 1b as d — > 0, we have 

lim sup II^ (ui , vi, u 2 , V2) 2 < 7r(supp7r + (ui — w 2 , «i — ^2)) ■ 

n— >-oo 

The push-forward measure of the mapping (ui , Vi , U2, V2) i-> (ui — u-i , «i — V2) on M 2 is absolutely 
continuous with respect to A2, so our argument is complete if we show that 7r(supp7r — z) = 
for almost every z £ R 2 . But this follows from the assumption that A2(supp7r) = 0, since 

J 7r(supp7r — 2;)A 2 (d2;) = J A 2 (supp tt — z)n(dz) = J A 2 (supp7r)7r(d2;) = 0, 

where the first equality follows from Lemma 1.28 in [23] and the second from the translation 
invariance of the Lebesgue measure. □ 



4. Central limit theorem 



The proof of the central limit theorem, Theorem 2.12, is based on a chaos decomposition of the 
power variation, that is, representing it as a L 2 -convergent series of iterated Wiener integrals 
with respect to the white noise W. Then, we apply the limit theory for iterated Wiener integrals 
to establish convergence of finite-dimensional distributions. We will begin by recalling some key 
facts of the chaos decomposition and the related central limit theorem. 

4.1. Central limit theorem via chaos decompositions 

Let us denote by T-L the Hilbert space L 2 ([—l, l] 2 ), which will have a special role in what follows. 
Moreover, let H® k = L 2 ([—l, l] 2fe ) be the fc-fold tensor product of H, for k 6 N, and denote by 
H &k the set of symmetric functions belonging to W® , that is, for any / g "H 0fc , permutation 
s : {1, . . . , k} — > {1, . . . , fc}, and almost any (zi, . . . , z^) £ [0, l] 2k , 

/(«!,. . . ,Zi) = /(2 s (l), ■ ■ ■ , 2 s (fc))- 

For any / £ W. &k , the fc-fold iterated Wiener integral of the kernel / with respect to the white 
noise W, denoted by Ik(f), can be defined as a linear map H &k — > L 2 (ilw) with the key 
property 

Vw[h(f) 2 ] =kl\\f\\ 2 um . 

(For the details of the construction, see [28, pp. 7-10].) The remarkable feature of these integrals 
is that any X £ L 2 (Qw) admits a unique chaos decomposition [28, Theorem 1.1.2], 

00 

X = J2 I k(fk) mL 2 (n w ), (4.1) 

fc=0 

where fk £ H &k for any fc £ No, with the convention that fo := Ew[X] and Iq is the identity 
map on R. 

If we are given a sequence of random variables in L 2 {Q,\y) and we would like to show that 
they converge in law to a Gaussian distribution, the chaos decomposition (4.1) turns out to be 
instrumental. Specifically, such convergence can be established by verifying some straightforward 
criteria on the associated kernel functions. To formulate the criteria, recall that for any r £ 
{1, . . . ,fc-l}, the r-th contraction of f<U = /f 5 ®-- ■®fj: 1) £ H m and = f[ 2) ®-- ■<®f < j? ) £ 
H® k is the function 

®r / (2) - II </£U fFXfP ® - fi-r » /ffi ® • • • ® /f e 

i=l 
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The following multivariate central limit theorem is a slight reformulation of Theorem 1 in [8], 
originally a corollary of the results of Nualart and Peccati [29] , and Peccati and Tudor [30] . 

Lemma 4.2 (CLT via chaos decompositions). Let d £ N and for any n E N, let x[ n) , . . . , X^ e 
L 2 (rt\y) be such that for any i = 1, . . . , d, 

oo 

fc=i 

where E H &k . Suppose that 
(1) for anyi = l,...,d, 

OO 

hn^limsup^ k\\\f$\\ 2 k = 0, 

III fU<J n ^QQ 

fc=m 

£/iere exist positive semidefinite d x d-matrices '^,,'^S l \^ 2 \ . . . such that for any i, j = 
l,...,d and k E N, 

and that J2T=i ^ = N , 
(3) for any % = 1, . . . , d, fc £ N, and r = 1, . . . , k — 1, 

^ \\f$ ®r fk"i\\n®HH-r) = 0. 

Then, (x[ n \. . .,X ( d n) ) 4 N d (0, N) as n -)• oo. 
4.2. Proof of Theorem 2.12 

Throughout this section, apart from formula (4.19) below, we work conditional on the realization 
of a, regarding it as deterministic — similarly to the earlier proof of Theorem 2.6. 

We introduce some convenient notation. We define for any n E N, and i, j = 1, . . . , L^ 1 ]) 
function f n ,(i,j) S "H by 

/n,(i,i)( s >*) : = h n {e n i - s,e n j - t)a^ s t ^ 

and its normalized counterpart f n ,(i,j) : = Wfn^i^W^ fn,(i,j) G H. By definition, Y (^ nil k n j)) = 
h(fn,(i,j))- Thus, we have 

^[Y(4£iM))*\=Mn,M)\& 

and 

Corr w [r(i?£ 2i , fcwi) ),r(i?£ l2 >knj2) )] = </n ) (M, il ),/n ) ( i2 ,i 2 )> w , (4-3) 

whence 

|(/n,(ii /n,(«2j2))«| (4-4) 

We also write 

Z U) ■■= e n c-^ 2 (V^ t) (k n ,n) - E w [V ( %(k n ,n)}), (s,t) E [0, l] 2 . 

As a preparation for the proof of Theorem 2.12, we prove some key lemmas. First, we obtain 
a uniform estimate for the decay of correlations (4.3) under Assumption 2.8. 
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Lemma 4.5 (Correlation estimate). If Assumption 2.8 holds, then 
Proof. For any (ix,ji) 7^ (« 2j 32), we have the bound 



p n := SUp K/n,(ii,Ji))/n,(» 2t jj))^| ^ °( £ n)- 

(il,il)^(»2,j2) 



(/n,(ii,ii): /n,(i2j2))^ 



< CT (y^(e,T) 1/2 7 r n (e + e„( i i-i 2 ),T + e„(j 1 -j 2 )) 1/2 dedr^ . (4.6) 
Since £7„ fl (-E„ + (e n i,e n jj) is a A 2 -null set for any (i, j) eZ 2 \ {0}, we have 



/(z)A 2 (d*)< / /(z)A 2 (dz) + / /(z)A 2 (dz). (4.7) 

R 2 R 2 \-E„ R 2 \(B„ + (e„(»2-»i),e„02-Ji))) 

for any / € L 1 ^ 2 ). By applying (4.7), using the inequality (s + t) 2 sC 2(s 2 + t 2 ), (s,t) G M 2 , 
the Cauchy-Schwarz inequality, and making the obvious change of variables, we find that the 
right-hand side of (4.6) is bounded by 

4TT n (R 2 \E n ) = o(e 2 ), 

which is independent of (ii,ji) and (i 2 ,j 2 ). □ 
Next, we derive a chaos decomposition for the random field Z^ n '. 

Lemma 4.8 (Chaos decomposition). For any n € N and (s,t) e [0, l] 2 , 

00 

41)=£ J *(*fe?) mL2 ^w), (4.9) 

k=2 

where 

[s/e n ] [t/e„i . || , ,|2 v p/2 

p (n,fe) ._ Ofe V"^ V"^ f \\}n.{i,j)\\u \ 7®k c o/igife 

— fe! £ « 2^ 2^ I c J fc rt • 

i=i 3=1 \ " / 

Proof. Since j fe 3 )) ~ N(0,l) given cr, and since the Hermite rank of the 

function u p is 2, we have the expansion 

41)=^ P/2 E E (l^«a„i))r- E ^[l y «a„i))n) 
1=1 J=l 



^ p/2 E E n/«,(«)iiw(ii/»,(ij)ii« 1 i r K£i, fcl j))) ( 41 °) 

i=l 3=1 

^- p/2 E E ii/«,fe)ii«E^(ii/».(^)ii« r Kw 



As H^I^i?^ \ k j)) = ^(fn^ij)) i the Hermite representation of iterated Wiener inte- 

grals [23, Theorem 13.25] yields 

By plugging (4.11) into (4.10) and rearranging, we arrive at the asserted chaos decomposition. 

□ 
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Remark 4.12. Since CC2,(X3, ■ ■ • are the non-zero coefficients in the Hermite expansion of u p , we 
have 

X) IT = / M x ) 2 "f( dx ) = ™2 P - rrip < oo. (4.13) 

k=2 

We will use Lemma 4.2 to prove the convergence of the finite-dimensional distributions of 
using the chaos decomposition (4.9). To this end, we study the asymptotic behavior of 
the kernels in (4.9). 

Lemma 4.14 (Asymptotics of kernels). If Assumption 2.8 holds, then for any (si,ti) , (^2^2) € 
[0, l] 2 , fc > 2, and r = 1, . . . ,k - 1, 

OO 

lim hmsup J2 ^\\ F l:;%\\ 2 n^ = °. (4.15) 

S1AS2 — so tiA*2 — to 

lim =^ / / a 2p , d M d«, (4.16) 

— so —*o 

J™ H^'S) ®r ^&,S)llUc^ = 0. (4.17) 

Proof. Below, we use the index sets 

Z„ := {£1 : 1 < ii < [si/e„\} x {ji : 1 < ji < Ui/ £ «J} x {i 2 : 1 £2 < |_ s 2/ £ nJ} 
x {j 2 : 1 < j 2 < \t 2 /e n \} 

and 

X n :=Z* \ L(«x A s 2 )/e n J, 1 < j < L(*i At 2 )/e rl J}. 

Let us expand 

L\/jp(n,k) p(n.k) \ 
\ (si,i 2 )' (S2:t 2 )/«® fc 

,(iiiji)llw HA, (»2,j2)im \ / \ 

I ^2 ) Vn,(n ji)' Jn,{iz,h)/H® h > 

s. a-\<zt V « / 



(u >ii,i2 ,j2)ex, 

wh ere </^ iliJ0 , % A) ) W .* = </n,(^i), /n.fe „•,)>*. As k > 2, we have by (4.4), 

SU P \{fn,(i 1 ,ji)i fn,(i 2 ,j2))-}{\ ~ SU P | (/", (il ,ji ) ' fn,(i 2 J2) | 

(»l,Jl>j2,32)eii (il ,31,12 j2)ex„ 



< SUp |(/n,(» 1 ,ji))/n,(i 2 ,i2)) 

(ii,3iiW,j2)eX 
2 / 2 



« 1 



« = o(<). 

Due to the bound |X n | ^$(si,ti),(s2,t2) e n 4 ' we can wr ^ e 

2 / L(siAs 2 )/EnJ L(*iA«2)/enJ ,u f \\2 \P \ 
w/p(«.' s ) pi™.' 8 ) \ _ a k\Jl \ " \ ^ / IUn,(i,j)ll^ \ p. \ 
fc! \^i,t2)'^( S 2,t2)/«^ - fcfl £ « <L 2s I " J + B fc,n|> 
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where lim n _ J . 00 sup fc ^ 2 |©fe,n| = 0. Now (4.16) follows since 

L(siAs 2 )/e„J L(*lAt 2 )/e„J 



el 



E E 

i=i j=i 



/ II 11?^ \ P 



LOiAs 2 )/ E „J L(*iAt 2 )/ £ „J , . j 

£ ™ E E ( /4«i-e,e-i-r) 7r n( d f' d ' r )) 

i=l 3=1 W 



siAs 2 ftiAta 

cr 2p ^dud?;, 



(cf. the proof of Theorem 2.6). Moreover, by the reverse Fatou's lemma and (4.13), we obtain 



hmsup £ < E t // // VV Mdw 0. 



JO 



establishing (4.15). 

To show (4.17), we may use the bound 



(«i,ti) ^ 1 (ai.tOllwaaCfc-r) 

(ii,i 2 ,i3,i4)ej; 



< 4 / f ® fe <X, / r<8lfe f®' 6?) f® fe \ 

~cr,/c t„ Vn,ii ^ r • / ri,i 2 ' •'"43 Wr •'n,i4/-H®2(fc-r) 



— ^ ' i /n,i 2 )^(/n,i3 > /«44 )-^(/n,ii ; /n,i3 )^ (fn,i 2 > /n,i4 i 

(ii,i 2 ,i3,i4)£Jn 

where 

J-„ := ({1, . . . , Lsi/e„J} x {1, ... , LtiAnJ }) 4 - 
Since r > 1 and fc — r ^ 1, we have by (4.4), 

&n ^ ^ (/n,ii ; /n,i 2 (fn,i3 7 In,i4 )^ (/n.ii ) fn.13 (/n,i 2 7 fn.14 

(ii,i 2 ,i3,i 4 )6 

^ £ n ^ ] I (/n,ii ) /n,ia)^ 1 1 (/n,i3, fn,ii)ji \ | (/n,ii ) fn,is)^i \ | (/n,i2J Jn,ii)n | (4-18) 

(ii,i2,i3>i4)e 

We use a simple combinatorial argument to deduce that the right-hand side (r.h.s.) of the 
inequality (4.18) tends to zero. To this end, we define M n : J7„ — ► {0, 1, 2, 4} by 

«n(il, 12, 13,14) = l{i 1= i 2 } + l{i 3 =i 4 } + l{ii=i 3 } + 1 {i 2 =i4>: 

where the value 3 is, indeed, never attained. It is straightforward to check that 

< (si)tl) £ - 8 , i^'ai})! <(. lltl ) e- 6 , 

KH{2})\ <( Sl , tl ) e; 4 , K 1 ({4})| < (si , tl) e-\ 

for all n 6 N. Hence, splitting the summation on the r.h.s. of (4.18) by 

E=E + E + E + E 

^ ^({O}) ^({l}) Un\{2}) «^ 1 ({4}) 

and using Lemma 4.5 leads to the bound 

r.h.s. of (4.18) < ff , fe , (si , tl ) 4( £ - 8 o(4)+£n 6 o(4)+^V4)+£n 2 ) = o(l). □ 
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We are now ready proceed to the actual proof of Theorem 2.12, building on the preceding 
three lemmas. 

Proof of Theorem 2.12. In this proof, unlike in the rest of the paper, the space D([0, l] 2 ) is 
endowed with the Skorohod topology (see Appendix B). Ultimately, we switch to the uniform 
topology using Skorohod coupling [23, Theorem 4.30] and Lemma B.3, as S( p ) is a continuous 
random field. 

Step 1: Reductions. It is clearly sufficient to show that 

z(n) gfe) injD (r 0) i]2) 5 

n— >oo 

where S P = (m 2p — m 2 ) 1 / 2 ^). The quadrant Brownian sheets (W^t)) (x t)e[o i] 2 ' ^ or * = 
1,2,3,4, defined by 

w (sl) - h(l[o,.Mo,t]), Wg } t) := /i(l[- s , ]x[o,t]), 

^2) : = J i(l[-.,(JIx[-t > o]) s ^(2) : = ^l(l[0, S ]x[-t,0]), 

(modulo taking continuous modifications) generate the c-algebra J^w- Thus, by Lemma C.l, it 
is sufficient to show that for any continuous, bounded test function ip : C([0, l] 2 ) 4 x C([— 1, l] 2 ) x 
D([0, l] 2 ) — ► K, 

E[^(ffW, W^, W< 3 >, WW, a, Z<»>)] — E[p(wW S°°)] , (4.19) 

where cr is stochastic. But, in the view of Fubini's theorem, it is clear that (4.19) follows if we 
simply show that 

(WW,W<V,WW,WW,ZW) — ^ (W^,W^,W^,W^,E ip) ) (4.20) 

n— >oo 

in C([0, l] 2 ) 4 x Z?([0, l] 2 ), with the realization of a kept fixed. As usual, we will prove (4.20) by 
establishing convergence of the finite-dimensional distributions first, and then showing tightness. 

Step 2: Convergence of finite- dimensional distributions. We fix arbitrary d € N and (s,t) := 
((si,ti), . . . , (sd,td)) € [0, l] 2d . Let us denote for any i = 1,2,3,4, 

(s,t) V (si,*i)' ' (s d ,td)l 

and for any n € N, 

7 ( n ) /T-(n) y(n) \ 

*(a,t) — \ Zj (s 1 ,t 1 )>---' Zj (s d ,t d )J- 



We would like to show that 



V (s,t)' (s,t)> (*,*)' (s,t)' (a,t)/ oa \ 



* (2) 



(4.21) 



where is the covariance matrix of (W^ t y wf®^, W^ t s,W^\^) and 



(s,t)' (s,t)' (s,t)' (s 



-to 

Note that w/ 1 ^, w/ 2 L, w/ 3 L, and are vectors of first-order Wiener integrals, whereas 



(n) 

the chaos decompositions of the components of Z^ g t ^ do not have any contributions from first- 
order integrals. Thus, (4.21) holds whenever the chaos decompositions of z9f t y n 6 N satisfy 
the conditions of Lemma 4.2, which is indeed the case due to Remark 4.12 and Lemma 4.14. 
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Step 3: Tightness. Since marginal tightness implies joint tightness, it suffices to show that 
the sequence (Z' n ') ngN is tight in D([0, l] 2 ). Let us write 

T n ■= {e n i : i = 0, 1, . . . , \_e~ 1 \ } U {1}, n £ N. 

Moreover, let R be a rectangle with vertices in T 2 , that is, R = (sijSa] x C^ijfe] for some 
si, s 2 ,ti 1 t 2 £ 7^ such that si < S2 and ii < £2- By [15, pp. 1658, 1665], it is sufficient to show 
that 

^ w [Z^\Rf] <<r, p \ 2 (R) 2 . (4.22) 

Recall that 



Z^(R) = z\ n \,-Z\ n \ . 

V 7 (S2,t 2 ) («2,*l) 

LS2/£„J 

p/2 



, 7(") 

LW^J 



E E H/n.(iJ)lft«p(ll/»,(«)ll« 1 l r (^£ 



i=L*i/e n J + lJ = L*i/£"J+ 1 



Since p n — > 0, there exists no G N such that sup„^„ p n < 1/12. Thus, by Lemma 4.23, below, 
we have for all n J? n , 



E W [Z^(R) 4 ] < a 4E 



|t a /e»J 



i=Lsi/enJ+li=L*i/e»J+l 

< p eUe- 8 \2(R) 4 p 4 n +e- 6 \ 2 (R) 3 p 2 n + s- 4 \ 2 (R) 2 ) 
^\ 2 (R) 2 (s- 4 p 4 n + e- 2 pl + l), 

where we used (3.10) and the inequalities {[s 2 /e n \ — [si/e n \)([t 2 /e n \ — [ti/e n \) ^ e~ 2 A2(i?) 
and X 2 (R) ^ 1. Finally, ~p n = o(e n ) implies that 



sup (e n 4 pi + e n 2 p 2 n + l) < 00 



whence (4.22) holds. 



□ 



It remains to prove the following moment estimate that we used above to establish tightness. 
It is similar to — albeit much less general than — Proposition 4.2 of [37]. The key difference, 
however, is that unlike in [37], here the underlying Gaussian random variables do not form 
a stationary process. (Alas, the assumption of stationarity renders Proposition 4.2 of [37] 
unapplicable in our setting.) The proof relics on a product moment bound due to Soulier [3G, 
Corollary 2.1] and a simple combinatorial argument. 

Lemma 4.23 (Fourth moment). Let [Xx,...,X n ) be a Gaussian random vector such that 
E[Xi] = and E[X 2 ] = 1 for all i = 1, . . . ,n. If f £ L 2 (R,j) has Hermite rank r £ N and 
p := sup iytj |E[XjX,]| p* for some p* £ (0, 1/12), then 



E 



E/( x *) 



< 



Proof. We use first the trivial bound 

1-1 



E 



(4.24) 
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By Corollary 2.1 of [36], we have 

\nf(X n )f(X l2 )f(X l3 )f(X u )}\ < f<p . p™n(n,i 2 ,i 3) u)/2 ) (4 25) 

where 

u n {ix,i 2 ,i3,ii) : = \{i '■ i = ik for some k and i ^ ii for any i ^ fc}|. 

is the number of unrepeated indices in i 2l i%, ii) (cf. the function u n in the proof of Lemma 
4.14). Note that u n is a mapping from {1, . . . , n} 4 onto {0, 1, 2, 4}, since it is impossible to have 
exactly three indices that are not repeated. It is key to observe that 

lu-HmiZn' + nZn 2 , fc 1 ({1})| < n 2 , 

fc X ({2})|<n 3 , S'PKn', 

for all n € N. (In the CcLS6 1/"^ (ii, i 2 , 13,44) = 0, either i\ = i 2 = h = 14 or there are two distinct 
pairs of repeated indices.) Thus, by (4.24) and (4.25), we obtain 



E 



J2m 



< f>p . n 4 p 2r + n 3 p r + n 2 p r ' 2 + n 2 < n 4 p 2r + n 3 p r + 



n 2 , 



which completes the proof. □ 

5. Asymptotics of concentration measures 

In this section, we prove Proposition 2.15 by deriving polynomial estimates for the integrals of 
the weight function g over some certain decisive subsets of M 2 . As the proof of Proposition 2.18 
uses a closely related argument, we merely sketch its main points. 

5.1. Proof of Proposition 2.15 

Let g € L 2 (M. 2 ) be given by (2.14). It will be convenient to consider the non-normalized measures 

[i n (ds,dt) '■= h n (s,t) 2 dsdt, n G N. 

Note that the support of p, n , like 7r„, is contained in [0, 1 + 1/n] 2 . Under the present assumptions 
g is a symmetric function, which clearly implies that also h n is symmetric. Thus, as p n is 
absolutely continuous with respect to the Lebesgue measure, we have 

, ro2v ^ ^((0,l + l/n) 2 \(0,£ ra ) 2 ) p n (T n \(0,e n ) 2 ) 

TT n {K \ b n ) - j— - — — , (5.1) 

£t n ((0, 1 + 1/n) 2 ) Pn(T n ) 

where T n := {(s,t) :0<i<s<l + 1/n}. On certain subsets of T n , the expression for 
h n (s,t) can be simplified significantly. To make use of this fact, we define for any n e N, 
E n ■= {(s,t) : < t < s < 1/n} and 

S r W : = ( e „,i) x (0,1/n), 

B n ] '■= {( s i0 : £ n < s < 1, s - 1/n < t < s}, 

:={(s,t) :e n <s< 1 + 1/n, l/n<t< (s-l/n)}, 
B {4) := (1,1 + 1/n) x (0,1/ti) U {(s,t) : 1 < s < 1 + 1/n, s-l/n <t< a}. 

(See Figure 1.) To prepare for the proof of Proposition 2.15, we establish next polynomial 
estimates for the asymptotic behavior of the /i n -measures of some decisive subsets of T n as 
n — > 00. 
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1+i 



/b [a) 



'o(2) ,/ 



B 



(3) 



o(l) 



1 £n 1 1 , I 

n n 

Figure 1: The sets £?„, B< 1} , B { n \ B ( n } , and si 4) . 



Lemma 5.2 (Polynomial bounds). Suppose that the assumptions of Proposition 2.15 hold and 
denote 

n := inf{n > 4 : k n ^ 2, |£(a;)| > /or aZZ x G (0, 1/n)}. 

Then, 



(1) 


Mn((0,£r 


) 2 nr„) > Q/ n-^ 1 - 


(2) 






(3) 




= 0(n- 3+K< - 2a +V), 


(4) 


Mn(si 3) ) 


= for all n ^ uq, 


(5) 


M„K 4) ) 


= o(n- 2 ). 



Proof. Throughout the proof, we denote /(s) := s a £(s) for all s € (0, 1). It is straightforward 
to check that 

'/(*), (s,t)eE n , 
/(a) -/(a -1/n), (M)6B» W , 



h n (s,t) = < 



/(«- 1/n) -/(<), ( s ,t)ei?i 2) , 
0, (s,t)esi : 



(5.3) 



(3) 



(1) The inclusion C (0,£ n ) 2 n T„ implies that /J, n (E n ) ^n((0,£«) 2 n T n ). By (5.3), we 



have 



where /(s) 2 ^ s~ 2a mf u€{0il/n) £{u) 2 > e s~ 2a . Thus, 



1/n ps rl/ n 

/ dt/( S ) 2 ds = / S /(s) 2 d S , 
Jo Jo 



1/n flf 71 

Sf(s) 2 ds> e / 8- 

o Jo 



l ds = 



u -2(l-a) 

2(1 - a) 
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(2) Due to (5.3), we may write 



"l/n pi 

H n {B^)= I dt / (/( S )-/(.s-l/n)) 2 d S 

(5.4) 

(f(s)-f(s-l/n)) 2 d S . 



(i 



By the mean value theorem, for any s £ (e„, 1), there exist £ s € [s — l/n, s] such that 

f(s) - f(s - l/n) = i/'fe), (5.5) 

where 

/'(^) = -<- a - 1 ^ s ) + C^'(6)- 

Thus, we have the bounds 

f(Zsf^2(a 2 C 2{a+1) \ml+€ 2a ¥'\\lo) 

ZaAl+ZsKs 2 ^ (5-6) 

By plugging (5.5) into (5.4) and applying (5.6) we arrive at 

i r°° 

- {f{ S )-f{s-l/n)) 2 &s< aI n-* s- 2 ^ds 

n Je n ie„-l/n (5.7) 

(3) Proceeding as above, we have by (5.3), 

vn(B^) = £ y'^ (f( S - i/n) f(t)) 2 dty s , 

where ^ t — (s — l/n) ^ l/n. Thus, 

(f(s - l/n) - f(t)) 2 < a . e n- 2 (s - l/n)- 2 ^ +1 \ 
by the mean value theorem and bounds analogous to (5.6). Moreover, 

/ (f( S -l/n)-f(t)) 2 dt d S < a ,,n" 3 / S ~ 2 ^ds, 

r„ \Js-l/n J Je n -l/n 

and the assertion follows from (5.7). 

(4) Obvious, by (5.3). 

(5) The estimate follows by observing that 

u n (B^) < X 2 (B^) sup MM) 2 , 
( S ,f)efli 4) 

where ^(B^) ^ 2/n 2 and sup^ s t ^ eB w h n (s,t) 2 — > as n — > oo because of the boundary 
condition lim s _j.i_ i(s) = 0. □ 
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Proof of Proposition 2.15. As mentioned above, Assumption 2.8 (in any form) implies Assump- 
tion 2.5 with 7r = S Zo . Thus, in view of (5.1), it suffices to show that 



£ U \ E n) = \ , n 7^ TTn r-.rr, \ > : 



e n 2 fi n (T n \ (0,£») 
Mn (T„ \ (0, e„) 2 ) + /i n ((0, e„) 2 n T n ) n^oo 

which is equivalent to 



At„(T„ \ (0,£„) 2 ) 



By Lemma 5.2 and the assumption s n X n K , 

4 

//„(T„\(0,e„) 2 ) ^5>„(£# 



1=1 



_ fo(n- 3+K ( 2Q+1 )), K G [l/(2a + 1), 1), 
[o(n- 2 ), KG (0,l/(2a + l)), 

and 

e^((0,e„) 2 nT n ) > a , t n- 2(1+K - a \ n > n . 
Thus, when k E (0, l/(2a + 1)) we have 

n 2(a-«) 

M„ > Q>/ — , 
o(l) 

whence (5.8) holds if n ^ a. In the case k € [l/(2a + 1), l), 

n 2(o-«-l) 



M n > ai< 



0( n -3+K(2a+l)) 



and, consequently, (5.8) holds provided that k < (2a + l)/(2a + 3). It remains to note that for 
a 6(0,1/2), 

2a + 1 1 , , 

a < s ^ < ^ 7, (5.9) 

2a + 3 2a + 1 v 7 

whereas for a € [1/2, 1), 

1 2a + 1 . 

s$ — a. (5.10 

2a + 1 2a + 3 y J 

The sufficiency of the asserted conditions can now be verified using the inequalities (5.9) and 
(5.10). □ 

5.2. Sketch of the proof of Proposition 2.18 

Let now g G L 2 {R 2 ) be given by (2.17). We redefine the sets E n , B { n\ , and B { n\ 

n £ N as indicated in Figure 2. Moreover, the measures fi n , n € N, are redefined accordingly. 

Proof of Proposition 2.18 (sketch). Similarly to the proof of Lemma (5.2), we have 

rl/n nljn 2(l-a) 
»n{E n )= tf(t) 2 dt> £ t - 2a+1 dt=— 

Jo Jo 2 (! - a ) 

and 

Vn{B {2) ) = - f {f{t) f{t l/n)) 2 At = 0{n-W*+*>). 
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1 + i 



>*A 



! :/ /' /' BO) 

i i i r> n 



1 — f 



11 + 1 

2 2 T ii 



H h 

1 1 + 



Figure 2: The redefined sets £„, B^, b£ 2) , B ( n 3) , and si 4) . The sets labeled with zeros are 
fi n -rm\\ sets for large n under the specification of g by (2.17). 

Additionally, ^„(£?i 4 '') = o(n~ 2 ). For the remaining two sets, we obtain 

Mr, 



" ' n Je n /2 



72 



and 



1 



1 



f(t- l/n) 2 dt <i - / r 2a &t, 

n Je n /2 n Je n /2-l/n 



and observing that £„/2x£„/2-l/nxe„xii K (note that x is an equivalence relation), we 
have 

Mn (B«) +Ai „(Bi 3 )) = O^-^^" 1 )). (5.11) 
To prove the assertion, it suffices to show that 

M n := ^n(^n) ( 

Since the contribution of (5.11) is dominant in the denominator of (5.12), we have 

M n > i<a nCta-D-^+i), 
whence (5.12) holds provided that n < (2a - l)/(2a + 1). □ 
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A. Uniform convergence of functions of two variables 



It is a well-known fact that if non-decreasing functions on [0, 1] converge pointwise to a contin- 
uous function, then the convergence is, in fact, uniform. In the proof of Theorem 2.6 we invoke 
the following analogous result that applies to functions on [0, l] 2 . 

Lemma A.l (Uniform convergence). Let /i, /a, • • • be functions [0, l] 2 — > K such that for any 
n e N. 

fn(s,t) < f n (u,v) ifsi^uandt^v, 
and let f G C([0, l] 2 ). If f n (s,t) — > f(s,t) for any (s,t) € [0, l] 2 ,, then f n -> / uniformly. 

Proof. The assertion follows from a straightforward adaptation of the standard argument used 
in the univariate case (see, e.g., [17, pp. 113-114]). □ 



B. On the two-variable generalization of the cadlag property 

We will review briefly the natural generalization of the cadlag property (continuity from the 
right with finite limits from the left) for functions on [0, l] 2 , following the formulation due to 
Neuhaus [27]. To this end, we introduce for any (s,t) £ [0, l] 2 the quadrants 

Qi(«,t) : = (5,1] x (t,l], Q 2 (M) := [0,s) x (t,l], 

Qs(s, t) := [0, s) x [0, t), Q A {s, t) := (s, 1] x [0, t). 

The space D([0, l] 2 ) consists of functions / : [0, l] 2 — > K such that for any (s,t) € [0, l] 2 , the 
following two conditions hold. 

• We have f(s n ,t n ) — > f{s,t) if (s n ,t n ) is a sequence in Qi(s,t) such that (s n ,t n ) — > (s,i), 

• For any i = 2,3,4, there exists fi(s,t) € K that satisfies f(s n ,t n ) —> fi(s,t) if (s„,t„) is a 
sequence in Qi{s,t) such that (s n ,t n ) — ¥ (s,t). 

In other words, / G D([0, l] 2 ) is continuous from the direction of the first quadrant (cf. the cad 
property) and has limits from the directions of the other three quadrants (cf. the lag property) . 
Clearly, we have C([0, l] 2 ) C D([0, l] 2 ). 

The space D([0, l] 2 ) can be endowed with the generalized Skorohod topology defined by 
Bickel and Wichura [15] and Neuhaus [27], which can be characterized in terms of convergence 
of sequences as follows. Let us denote by A the class of mappings A : [0, l] 2 — > [0, l] 2 such that 
X(s,t) = (A (1) (s), A (2) (i)), where A (1) and A (2) are increasing bijections [0, 1] — > [0, 1]. 

Definition B.l (Skorohod topology). Let /, ft, / 2 . . . £ D([0, l] 2 ). We say that /„ — > / in the 
Skorohod topology if there exist Ai, A2, • • • € A such that 

sup |/„oA„(s,t)-/M)l+ sup ||A„(s, t) — (s, t)|| >0. (B.2) 

( s ,i)e[o,i] 2 ( s ,*)e[o,i] 2 ™^°° 

There is a Skorohod metric on D([0, l] 2 ) that is consistent with the convergence defined 
above (see [27, p. 1289]). Equipped with this metric, D([0, l] 2 ) enjoys the usual properties of 
separability and completeness (i.e., it is a Polish space), similarly to D([0, 1]). 

We use the Skorohod topology merely as a technical tool to establish convergence in law 
in the proof of Theorem 2.12, using the relatively tractable tightness criteria for the Skorohod 
topology [15, pp. 1665-1666]. Since the limit obtained in Theorem 2.12 is a continuous random 
field, we may — equivalently — equip D([0, l] 2 ) with the (non-separable) uniform topology, 
thanks to the following result. 
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Lemma B.3 (Uniform convergence). Let /i,/2,--. G D([0, l] 2 ) and / G C([0, l] 2 ). Tften, 
/„ — > / m £/ie Skorohod topology if and only if /„ — > / uniformly. 

Proof. It is obvious that uniform convergence implies convergence in the Skorohod topology. To 
show the converse, let us fix e > 0. Since / is uniformly continuous, there exists 5 > such that 
\f(s,t) - f(u,v)\ < e/2 if \\(s,t) - (u,v)\\ < S. Now, let Ai,A 2) ... G Abe such that (B.2) holds. 
Then there exists no G N such that for all n ^ no, 

sup \f n o Xn(s,t) - f{s,t)\ + sup \\X n (s,t) - (s,t)\\ < | A S. 
(s,t)e[o,i] 2 (s,t)e[o,i] 2 1 

By the triangle inequality, we have thus for all n ^ n , 

sup |/„(M) ~ /(»,*) I = sup |/„oA„(s,t)-/oA„(s,t)| 
( s ,t)e[o,i] 2 (M)e[o,i] 2 

< sup |/ n o A n (s,i) - /(s,t)| 
O,t)e[o,i] 2 

+ sup |/(s,t) -/o X n (s,t)\ < £, 
( s ,t)e[o,i] 2 

which completes the proof. □ 

C. Stable convergence lemma 

The following simple lemma is a key tool in proofs of stable convergence in law. It is certainly 
well-known and, indeed, used (implicitly) in several papers (e.g., [6, 8]), but due to lack of a 
reference, we provide a proof for the convenience of the reader. 

Lemma C.l (Stable convergence). Let U and V be Polish spaces. If U, U\, U2, ■ ■ ■ are random 
elements in U and V is a random element in V, all defined on a common probability space 
(Q', such that 

(U n , V) — ^ (U,V), 



then 



u n ^ u. 



Proof. We will use a monotone class argument. To this end, let / G C(U,M.) be bounded and 
write 

Mf.= \xe L°°(fi', F',P') : lim E'[f(U n )X] = E'[f(U)X}\. 

Clearly, ftAf is vector space that contains all constant random variables. Moreover, if X G 
L°°(il, J",P) and X G F, then 

\E'[f(U n )X]-Ey(U)X}\ < f \E'[f(U n )X] - Ey(U)X]\+E'[\X - X\}. 

Hence, M.f is closed under uniform convergence and if (X n ) C Aif is such that ^ X\ ^ X2 ^ 
■ • • ^ M for some constant M > 0, then linin^oo X n € Mf. Now, note that 

C := {ip(V) : ip G C(V,M) is bounded} 

is closed under multiplication and C C A4f by the continuous mapping theorem. Thus, by 
the functional monotone class lemma [20, p. 14], Mf contains any bounded o~(C) -measurable 
random variable. Since V is separable, we have a(V) = a(C) and the assertion follows. □ 
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